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found to be well described by four-dimensional Yukawa-type function e~~ mr / r in the infrared and 
intermediate region of r = [x^x^l 2 = 0.1 ^ 1 .Ofm. The infrared effective gluon mass is obtained 
asm~ 600MeV. Associated with the 4D Yukawa-type gluon propagator, we derive analytical 
expressions for the zero-spatial-momentum propagator D§{t ), the effective mass M e g(t), and the 
spectral function p (<b) of the gluon field. Remarkably, the obtained gluon spectral function p(co) 
is almost negative definite, except for a positive 5-functional peak at co = m. Since the Yukawa- 
type propagation indicates a three-dimensional space-time character, we consider a hypothesis 
of an effective dimensional reduction by generalized Parisi-Sourlas mechanism in a stochastic 
color-magnetic vacuum of infrared QCD. 
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1. Introduction 

The analysis of gluon properties is an important key point to clarify the nonperturbative aspects 
of QCD [[I], [2], ||] . In particular, the gluon propagator, i. e. , the two-point Green function is one of the 
most basic quantities in QCD, and has been investigated with much interests [[I], f|, |5j. Dynamical 
gluon-mass generation is also an important subject related to the infrared gluon propagation. While 
gluons are perturbatively massless, they are conjectured to acquire a large effective mass as the 
self-energy through their self-interaction in a nonperturbative manner [|,|]. Actually, glueballs, 
color-singlet bound states of gluons, are theoretically predicted to be fairly massive, e.g., about 
1.5GeV for the lowest ++ and about 2GeV for the lowest 2 ++ , in lattice QCD calculations 0. 

For the direct investigation of the gluon field, gauge fixing is to be done. Among gauges, 
the Landau gauge is one of the most popular gauges in QCD, and it keeps Lorentz covariance and 
global S\J(N C ) symmetry. In Euclidean QCD, the Landau gauge has a global definition to minimize 
the global quantity R = J d 4 x Tr{A^ (x)A^ (x)} = ^ J d A xA a ^ (x)A^ (x) by gauge transformation. The 
local condition d^A^ (x) = is derived from the minimization of R. The global quantity 7? can be 
regarded as "total amount of the gluon-field fluctuation" in Euclidean space-time. In the global def- 
inition, the Landau gauge has a clear physical interpretation that it maximally suppresses artificial 
gauge-field fluctuations relating to gauge degrees of freedom [jl|]. 

In lattice QCD, the Landau gauge is defined by the maximization of 7?i att = ReTrL^(x), 
with the link- variable U^(x) = e mgA ^ x ) {a: lattice spacing, g: QCD gauge coupling). The gluon 
field is defined as A^(x) = 2E^{Ufi(x) — uJl(x)} — (trace part). In the Landau gauge, the mini- 
mization of gluon-field fluctuations justifies the expansion by small lattice spacing a. In Euclidean 
metric, the gluon propagator is defined by the two-point function as D a ^ v (x — y) = {A a ^{x)A h v (y)). 
Here, owing to the symmetries and the transverse property, the color and Lorentz structure of the 
gluon propagator is uniquely determined in the Landau gauge. 

In this paper, using SU(3) lattice QCD Monte Carlo calculations, we study the functional 
form of the Landau-gauge gluon propagator, D(r) = ^ N \ ^ D a ^(x) = 3 ^ N l_^ (A^(x)A^(O)), as 
a function of 4D Euclidean distance r = (x^x^) 1 / 2 . We mainly deal with the coordinate-space 
propagator D(r) for the infrared and intermediate region of r = 0. 1 ~ l.Ofm, which is relevant for 
quark-hadron physics. Based on the obtained function form of the gluon propagator, we aim at a 
nonperturbative description of gluon properties, 

2. Functional form of Landau-gauge gluon propagator 

The SU(3) lattice QCD Monte Carlo calculations are performed at the quenched level using 
the standard plaquette action with j3 = 2N c /g 2 =5.7, 5.8, and 6.0, on the lattice size of 16 3 x 32, 
20 3 x 32, and 32 4 , respectively. The lattice spacing a is found to be a = 0.186,0.152, and 0.104fm, 
at p = 5.7, 5.8, and 6.0, respectively, when the scale is determined so as to reproduce the string 
tension as yfc = 427MeV from the static QQ potential |jg|. Here, we choose the renormalization 
scale at pL = 4GeV for j8 = 6.0, and make corresponding rescaling for j3=5.7 and 5.8 [JT]]. 

Figure 1(a) and (b) show the coordinate-space gluon propagator D(r) and the momentum- 
space gluon propagator D(p 2 ) = J d 4 x e ipx D{r), respectively. Our lattice QCD result of D{p 2 ) is 
consistent with that obtained in previous lattice studies, although recent huge- volume lattice studies 
[Q] indicate a suppression of the gluon propagator in the Deep-IR region (p < 0.5GeV). 
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Figure 1: (a) Lattice QCD results (symbols) of the Landau-gauge gluon propagator D(r) = D a ^ l [x)/2A in 
coordinate space, and 4D Yukawa-type function DyukawaM = Ame~'" r /r (solid line) with m = 0.624GeV 
and A = 0. 162. The dash-dotted line denotes a typical example of the massive-vector propagator D mass (r). 
(b) The Landau-gauge gluon propagator £>{p 2 ) in the momentum space of p^ — | sin(-^). The solid line 

denotes 4D Fourier-transformed Yukawa-type propagator, i.e., Syukawa(p 2 ) = 47Z 2 Am(p 2 + m 2 )~ 3 / 2 . 



We find that the lattice gluon propagator D(r) cannot be described by the free massive Eu- 
clidean propagator D mass (r) = / p^re _ip '*p^ = 4^7 K i ( mr ) CMz): modified Bessel function) 
[]|] in the whole region of r = 0. 1 ~ l.Ofm, as shown in Fig. 1(a). 

By the functional-form analysis, we find that the Landau-gauge gluon propagator D(r) in the 
coordinate space is well described by the 4D Yukawa-type function [[I]] 

1 e~ mr 
D(r)^-D a ^(r)=Am—, (2.1) 

with m = 0.624(8)GeV and A = 0.162(2) in the range of r = 0.1 ~ l.Ofm, as shown in Fig. 1(a). The 
gluon propagator D(p 2 ) in the momentum space is also well described by 4D Fourier-transformed 
Yukawa-type function as D{p 2 ) = ^D a ^(p 2 ) = (p2 4 ^ 3/2 for 0.5GeV < p < 3GeV [fj]. 

3. Analytical applications 

In this section, as applications of the Yukawa-type gluon propagator, we derive analytical 
expressions for the zero-spatial-momentum propagator D^{t), the effective mass M e ff(t), and the 
spectral function p(co) of the gluon field [jj]]. All the derivations can be analytically performed, 
starting from the Yukawa-type gluon propagator ^YukawaM- Although the real gluon propagator 
has some deviation from the Yukawa-type in UV region, this method is found to be workable to 
reproduce lattice QCD results, as shown below. 

3.1 Zero-spatial-momentum propagator of gluons 

First, we consider zero-momentum gluon propagator Do(t) = ^ Y.xi^i (-*^Mu (0> 0)) = L.v^( r ) > 
where r = Vx 2 + t 2 is the 4D Euclidean distance. For the simple argument, we here deal with the 
continuum formalism with infinite space-time. Starting from the Yukawa-type gluon propagator 
^Yukawa (r), we derive the zero-spatial-momentum propagator as |jl|] 

D (t) = [ d 3 x DyukwJV^ + t 2 ) = 47iAm f dx - % g-'"^^ 7 = AnAtK^ (mt). (3. 1) 
J Jo Vx 2 + 1 2 
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Figure 2: (a) The zero-spatial-momentum propagator Do (f ) of gluons in the Landau gauge, (b) The effective 
mass M e g(t) of gluons in the Landau gauge. The symbols are the lattice QCD data on 32 4 at j3 = 6.0, and the 
solid line is the theoretical curve derived from 4D Yukawa-type propagator with w=0.624GeV and A=0. 162. 



In Fig.2(a), we show the theoretical curve of DoO) in Eq.(|l]) with m=0.624GeV and A=0.162, 
together with the lattice QCD result of Dn(f) in the Landau gauge. For the actual comparison with 
the lattice data, we take account of the temporal periodicity The lattice QCD data are found to 
be well described by the theoretical curve, associated with the Yukawa-type gluon propagator. 



3.2 Effective mass of gluons 

Second, we investigate the effective mass M e ff(f) of gluons. The effective mass plot is often 
used for hadrons as a standard mass measurement in lattice QCD. For the simple notation, we use 
the lattice unit of a = 1 in this subsection. In the case of large temporal size, the effective mass is 
defined as M ef£ (t) = hx{D (t)/D (t + 1)}. 

In Fig.2(b), we show the lattice result of M e f[(t), where we take account of the temporal peri- 
odicity. The effective gluon mass exhibits a significant scale-dependence, and it takes a small value 
at short distances. Quantitatively, the effective gluon mass is estimated to be about 400 ~ 600MeV 
in the infrared region of about lfm [JXJ] _ This value seems consistent with the gluon mass suggested 
by Cornwall [Q], from a systematic analysis of nonperturbative QCD phenomena. 

Now, we consider the consequence of 4D Yukawa-type propagator DyukawaM of gluons. For 
simplicity, we here treat the three-dimensional space as a continuous infinite-volume space, while 
the temporal variable t is discrete. We obtain an analytical expression of the effective mass [Jl]], 

D Q (t) tKAmt) 
Meff(f) = ln — p-L^=\ n - _ (3.2) 

when the temporal periodicity can be neglected. In Fig.2(b), we add by the solid line the theoretical 



curve of M eS (t) in Eq.p^ with ra=0.624GeV. The lattice QCD data of M eff (f) are found to be 
well described by the theoretical curve derived from the Yukawa-type gluon propagator. From the 
asymptotic form K\(z) °< z~ l ^ 2 e~ z , the effective mass of gluons is approximated as §\§ 

Afeff(t)— »i-i]n(l + -) (for large t). (3.3) 
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This functional form indicates that M e ff(f) is an increasing function and approaches m from below, 
as t increases. Then, the mass parameter m ~ 600MeV in the Yukawa-type gluon propagator has a 
definite physical meaning of the effective gluon mass in the infrared region. 

Note that the simple analytical expression reproduces the anomalous "increasing behavior" of 
the effective mass M e ff(f) of gluons. Thus, this framework with the Yukawa-type gluon propagator 
gives an analytical and quantitative method, and is found to well reproduce lattice QCD results. 

3.3 Spectral function of gluons in the Landau gauge 

As a general argument, an increasing behavior of the effective mass M e n(t) means that the 
spectral function is not positive-definite More precisely, the increasing property of M e ff (?) can 
be realized, only when there is some suitable coexistence of positive- and negative- value regions in 
the spectral function p(co) However, the functional form of the spectral function of the gluon 
field is not yet known. 

The relation between the spectral function p{(o) and the zero-spatial-momentum propagator 
Do(?) is given by the Laplace transformation, Dq(/) = Jq dco p(ft)) e~ m '. When the spectral func- 
tion is given by a 8 -function such as p(ft)) ~ 5 (ft) — (Do), which corresponds to a single mass 
spectrum, one finds a familiar relation of D Q (t) ~ e~ m ". For the physical state, the spectral func- 
tion p((o) gives a probability factor, and is non-negative definite in the whole region of ft). This 
property is related to the unitarity of the S -matrix. 

From the analytical expression of the zero-spatial-momentum propagator Do(t) = AnAtK\ (mt), 
we can derive the spectral function p (ft)) of the gluon field, associated with the Yukawa-type gluon 
propagator [[I]]. For simplicity, we take continuum formalism with infinite space-time. Using the 
inverse Laplace transformation of the modified Bessel function, we derive the spectral function 
p (ft)) of the gluon field as [jl|] 

. . AnAm n . . AnA/\/2m . ,„ „ N 

P^ = -J772 2W2 e co-m-e ) + — S(<o-m-e), (3.4) 

(ft) z — m z yl z {(0 — m) i l A 

with an infinitesimal positive £, which is introduced for a regularization. Here, m ~ 600MeV is 
the mass parameter in the Yukawa-type function for the Landau-gauge gluon propagator. The first 
term expresses a negative continuum spectrum, and the second term a 5-functional peak with the 
residue including a positive infinite factor as £ -1 / 2 at ft) = m + £. 

We show in Fig.3 the spectral function p(ffl) of the gluon field. As a remarkable fact, the 
obtained gluon spectral function p(ft)) is negative-definite for all the region of ft) > m, except 
for the positive 5-functional peak at ft) = m. The negative property of the spectral function in 
coexistence with the positive peak leads to the anomalous "increasing behavior" of the effective 
mass M e ff(t) of gluons [jl|]. Actually, Eq^^l) leads to Eq.(3.2), which well describes the lattice 



result of the effective mass M e ff(f), as shown in Fig.2(b). 

We note that the gluon spectral function p(ft)) is divergent at ft) = m + e, and there are two 
divergence structures: a 5-functional peak with a positive infinite residue and a negative wider 
power-damping peak. On finite- volume lattices, these singularities are to be smeared, and p(ft)) 
is expected to take a finite value everywhere on ft). On the lattice, the spectral function p(ft)) 
is conjectured to include a narrow "positive-valued peak" stemming from the 5-function in the 
vicinity of ft) = m (+e) and a wider "negative-valued peak" near ft) ~ m in the region of ft) > m |jl|]. 
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Figure 3: The spectral function p (<b) of the gluon field, associated with the Yukawa-type propagator. The 
unit is normalized by the mass parameter m ~ 600MeV. p (to) shows anomalous behaviors: it has a positive 
5-functional peak with the residue of +°° at CO — m, and takes negative values for all the region of co > m. 

In this way, the Yukawa-type gluon propagator indicates an extremely anomalous spectral 
function of the gluon field in the Landau gauge. The obtained gluon spectral function p(co) is 
negative almost everywhere, and includes a complicated divergence structure near the "anomalous 
threshold", CO = m (+e). Thus, this framework with the Yukawa-type gluon propagator gives an 
analytical and concrete expression for the gluon spectral function p(©). 

4. Effective dimensional reduction in gluonic vacuum by Parisi-Sourlas mechanism 

We discuss the Yukawa-type gluon propagation and a possible dimensional reduction due to 
the stochastic behavior of the gluon field in the infrared region [Q]. As shown before, the Landau- 
gauge gluon propagator is well described by the Yukawa function in /owr-dimensional Euclidean 
space-time. However, the Yukawa function e~ mr / r is a natural form in ?/jree-dimensional Euclidean 
space-time, since it is obtained by the three-dimensional Fourier transformation of the ordinary 
massive propagator (p 2 + m 2 ) -1 . In fact, the Yukawa-type propagator has a "three-dimensional" 
property. In this sense, as an interesting possibility, we propose to interpret this Yukawa-type 
behavior of the gluon propagation as an "effective reduction of the space-time dimension". 

Such a "dimensional reduction" sometimes occurs in stochastic systems, as Parisi and Sourlas 
pointed out for the spin system in a random magnetic field |9p. In fact, on the infrared dominant 
diagrams, the D-dimensional system coupled to the Gaussian-random external field is equivalent 
to the (D — 2) -dimensional system without the external field, due to a hidden SUSY structure. 

We note that the gluon propagation in the QCD vacuum resembles the situation of the system 
coupled to the stochastic external field. Actually, as is indicated by a large positive value of the 
gluon condensate (G^ V G« v ) = 2(H 2 , — E 2 ) > in the Minkowski space, the QCD vacuum is filled 
with a strong color-magnetic field [JTOQ „ which can contribute spontaneous chiral-symmetry break- 



ing Ql lp, and the color-magnetic field is considered to be highly random at the infrared scale. Since 
gluons interact with each other, the propagating gluon is violently scattered by the other gluons in 
the randomly-oriented color-magnetic fields of the infrared QCD vacuum, as shown in Fig.4. 

Actually at the infrared scale, the gluon field shows a strong randomness due to the strong 
interaction, and this infrared strong randomness is considered to be responsible for color confine- 
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Figure 4: A schematic figure for a propagating gluon. The QCD vacuum is filled with color-magnetic fields 
which are stochastic at an infrared scale, and the gluon propagates in the random color-magnetic fields. 

ment, as is indicated in strong-coupling lattice QCD. Even after the removal of fake gauge degrees 
of freedom by gauge fixing, the gluon field exhibits a strong randomness accompanying a quite 
large fluctuation at the infrared scale. 

As a generalization of the Parisi-Sourlas mechanism, we conjecture that the infrared structure 
of a theory in the presence of quasi-random external fields in higher-dimensional space-time has 
a similarity to the theory without the external field in lower-dimensional space-time [JIJ]. From 
this point of view, the Yukawa-type propagation of gluons may indicate an "effective reduction of 
space-time dimension" by one, reflecting the interaction between the propagating gluon and the 
other gluons in randomly-oriented color-magnetic fields in the infrared QCD vacuum. 

In any case, it is an interesting and important subject to clarify the nonperturbative QCD 
vacuum structure in terms of gluonic properties [§] including the gluon propagation (Jxjl . 
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